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Abstract 

We construct an explicit 50(10) x U(1)f model of the Yukawa interactions by us- 
ing as a guide previous phenomenological results obtained from a bottom-up approach 
to quark and lepton mass matrices. The global U(l)p family symmetry group sets 
the textures for the Majorana and generic Dirac mass matrices by restricting the type 
and number of Higgs diagrams which can contribute to each matrix element, while the 
5O(10) group relates each particular element of the up, down, neutrino and charged 
lepton Dirac matrices. The Yukawa couplings and vacuum expectation values associ- 
ated with pairs of 1, 45, 10, and 126 Higgs representations successfully correlate all 
the quark and lepton masses and mixings in the scenario incorporating the nonadia- 
batic solar neutrino and atmospheric neutrino depletion effects. 
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In a series of manuscripts [1], the authors have demonstrated how a new bottom-up 
approach to the quark and lepton mass and mixing problem can be used to construct phe- 
nomenological quark and lepton mass matrices at the supersymmetric SO (10) grand unifi- 
cation scale, which lead to the assumed experimental input at the low scales. As such, this 
procedure provides an alternative to the usual top-down approach [2], where mass matrices 
are constructed based on some well-defined theoretical concepts. Of special interest is a set 
of mass matrices found by our approach which exhibit a particularly simple 5*0(10) structure 
for the scenario based on the depletions of solar [3] and atmospheric [4] neutrinos through 
oscillations. 

In this letter we construct an explicit SO (10) xU(1)f model at the grand unification scale 
by making use of the phenomeno logical mass matrices as a guide. The global U(1) F family 
symmetry singles out a rather simple set of tree diagrams which set the textures for the Dirac 
and Majorana mass matrices, while 50(10) relates the corresponding up, down, neutrino 
and charged lepton Dirac matrix elements to each other. The quantitative numerical results 
obtained from the model agree in detail with the input data assumed for the bottom-up 
approach. 

The starting point for our bottom-up approach was the reasonably well-known quark 
mass and Cabbibo-Kobayashi-Maskawa mixing matrix data [5]. To this we appended neu- 
trino mass and mixing data consistent with the nonadiabatic Mikheyev-Smirnov-Wolfenstein 
(MSW) [6] resonant matter oscillation depletion [3] of the solar electron-neutrino flux to- 
gether with atmospheric muon-neutrino depletion [4] through oscillations into tau-neutrinos. 
After running the Yukawa couplings to the grand unification scale, we applied Sylvester's the- 
orem, as illustrated by Kusenko [7] for quark data alone, to reconstruct complex-symmetric 
mass matrices. The construction is clearly not unique, but one can vary two parameters 
which determine the weak bases in order to select a set of mass matrices which exhibit par- 
ticularly simple 50(10) structure for as many matrix elements as possible. We refer the 
interested reader to Ref. [1] for details and begin here with the special phenomenological 
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matrices singled out by this procedure: 



M Nd 



M 



D 



M 



E 



diag(Y26; 126; 10) 
/ 10', 126 10', 126' 10' \ 
T26 10' 



(la) 
(16) 



10', 126 

V 10' 10' 10 J 

with M^, Mf 2 and M|[ anomalously small and only the 13 and 31 elements complex. Entries 

in the matrices stand for the Higgs representations contributing to those elements. Recall 

that the 50(10) product rules read 



16 x 16 = 10 s + 120 a + 126, 



16 x 16 = 1 + 45 + 210 



(2a) 
(26) 



We have assumed that vacuum expectation values (VEVs) develop only for the symmetric 
representations 10 and 126 and for 1 and 45. The Majorana neutrino mass matrix M R , 
determined from the seesaw formula [8] with use of M N ° irac and the reconstructed light 
neutrino mass matrix, exhibits a nearly geometrical structure given by [9] 

/ F -Vfe ^fFC \ 



M 



R 



-Vfe e 



-^Tec 



(3) 



V Vfc -Vec c J 

where E = \\/FC with all elements relatively real. It can not be purely geometrical, 
however, since the singular rank-1 matrix can not be inverted as required by the seesaw 
formula, M N "ff ~ -M^^^M^^M^rac. 

The challenge is now to introduce a family symmetry which will enable one to derive 
the mass matrix patterns in a simple fashion. For this purpose, we propose to use a global 
U(l)p family symmetry [10] and to reduce the problem to the construction of one generic 
Dirac matrix, M£n rac , along with the single Majorana matrix, M R . As noted above, the 
5*0(10) symmetry will relate the corresponding matrix elements of the four Dirac matrices 
to each other. The important roles played here by supersymmetry (SUSY) are twofold. Not 
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only does SUSY control the running of the Yukawa couplings between the SUSY GUT scale 
and the weak scale where it is assumed to be softly broken, but it also allows one to assume 
that only simple tree diagrammatic contributions to the mass matrices need be considered 
as a result of the nonrenormalization theorem applied to loop diagrams. This tree diagram 
procedure was first suggested by Dimopoulos [11] twelve years ago. 

Simplicity of the 5*0(10) structure requires that just one Higgs 10 representation con- 
tributes to the (Mr)i rac ) 33 element (hereafter labeled D33), i.e., we assume complete unifica- 
tion of the Yukawa couplings at the unification scale: fn T = m b = m t / tan(3 10 , where tan/3 10 
is equal to the ratio of the up quark to the down quark VEVs in the 10 

rh t = g w (v/ y/2) sin (3 W = g w v u 
rh b = m T = g w (v/y/2) cos(3 w = g w v d (4a) 
tan/3i =v u (5)/v d (5) 

in terms of the SU(5) decomposition of SO(10) with v = 246 GeV. The same 10 can not 
contribute to D23 = D32, for the diagonal nature of M u and M Nmrac requires the presence 
of another 10' with 

tan ft* = <(5')M(5') = (4b) 
Likewise we assume a pure 126 contribution to D22 with 

tanfe = w u (5)/w d (4E) (4c) 

The tree diagrams for these M Dirac matrix elements are illustrated in Fig. la. 

We shall now assign U{l)p charges to the three families (in order of appearance) and to 
the Higgs representations as follows with the numerical values to be determined later: 

16°, 16^, 16?, 10 a , I26 c (5a) 

Conservation of U (1) p charges then requires 2a + a = 0, a + (3 + b = and 2(3 + c = as 
seen from the diagrams in Fig. la. 
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In the above, we have taken the 2-3 sector of Mr, irac to be renormalizable with two 10's 
and one 126 developing low scale VEVs. We assume the rest of the M Dirac elements arise 
from non-renormalizable contributions with the leading ones shown in Fig. lb. For D13 we 
introduce a 45 e x Higgs field and construct an explicitly complex-symmetric contribution with 
the dimension-6 diagram, for which U(l)p charge conservation requires a + ^ + b + 2e = Q. 
This 45 x Higgs field develops a VEV in the direction which breaks 50(10) -> SU(5) x U(l) x 
with the SU(5) subgroup remaining unbroken. For D12 we introduce a different 45| Higgs 
field which breaks SO (10) — > flipped 577(5) x £7(1) and is related to the orthogonal 45 x and 
45y hyper charge VEVs by 

6 1 

< 45 z >= - < 45 x > ~- < 45y > (6) 

5 5 

as given in Table I. This Higgs field contributes to M® 2 but not to Mj|; thus it generates a 
zero in this position for the charged lepton mass matrix as suggested in (lb). The U(l)p 
charge conservation equation reads j3 + 7 — b + 2h = 0, as the 10'* Higgs field is required 
here to reduce the number of contributing diagrams. The Dll element is dimension-8 or 
higher and is left unspecified. The complex-symmetric Yukawa diagrams which we wish to 
generate are then neatly summarized by the ordering of the Higgs fields: 



£>33 : 


16 3 


-10 - 


16 3 






7J23 : 


16 2 


-10'- 


- 16 3 






7J32 : 


16 3 


-10'- 


- 16 2 






D22 : 


16 2 


-126 


- 16 2 






D13 : 


I61 


-45 x 


-10' 


-45x 


-I63 


D31 : 


16 3 


-45 x 


-10' 


-45 x 


-16! 


£>12 : 


I61 


-45 z 


-10'* 


-45 z 


-16 2 


£>21 : 


16 2 


-45 z 


-10'* 


-45 z 


-16x 



In order to obtain a different set of diagrams for the Majorana matrix, we begin the 



M33 contribution with a dimension-6 diagram shown in Fig. lc by including a new 126' 
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Higgs which develops a VEV at the GUT scale in the SU (5) singlet direction, along with a 
pair of l 9 Higgs fields. Here 2a + d + 2g = 0. The nearly geometric structure for M R can 
then be generated by appending more Higgs fields to each diagram. For M23 we introduce 
another 1'* Higgs field to construct a diagram with one 126^, one 45 x, one 1'^ and two l 9 
fields with charge conservation demanding a + (3 + d + 2g + e + f = 0. The new 1' field is 
needed in order to scale properly the Majorana matrix elements relative to each other. The 
remaining leading-order diagrams of the complex-symmetric Majorana mass matrix follow 
by appending more 45 e x , 45| and l' f Higgs lines. The pattern is made clear from the charge 
conservation equations: 2(3 + d + 2g + 2e + 2f = for M22, a + >y + d + 2g + e + h + 2f = 
for M13, (3 + 7 + d + 2g + 2e + h + 3f = for M12, and 2 7 + d + 2g + 2e + 2h + 4/ = for 
Mil. 

In summary, the following Higgs representations have been introduced in addition to 
those in (5a): 

I26 ,d , 453,, 45|, V, l' f (56) 

all of which generate massive VEVs near the GUT scale. In order to obtain CP-violation 
in the quark and lepton mixing matrices, we allow the VEVs for 45 x , 45^, 1 and 1' to be 
complex, but the VEVs associated with the 10, 10' , 126 and 126' representations can be 
taken to be real without loss of generality as seen from our bottom-up results. Clearly, many 
permutations of the Higgs fields are possible in the higher-order diagrams. 

At this point a computer search was carried out to generate U(l)p charge assignments 
leading to the fewest additional diagrams allowed by charge conservation. An especially 
interesting charge assignment stood out for which 

a = 9, (3 = -1, 7 = -8 

(8a) 

a = -18, b = -8, c = 2, d = -22, e = 3.5, / = 6.5, g = 2.0, h = 0.5 
One should note that since a + /3 + 'y — 0, the [S'O(IO)] 2 xU(1)f anomaly vanishes, whereas 
the [[/(l)^] 3 anomaly does not. Simplicity then suggests that the U(l)p family symmetry 
group be global with a familon being generated upon its breaking. 
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With the above charge assignments we can greatly limit the number of permutations 
and eliminate other unwanted diagrams by restricting the U(l)p charges appearing on the 
superheavy internal fermion lines. With the following minimum set of allowed charges for 
the left-handed superheavy fermions Fl and their mirror partners F£ 

F L : -0.5, 1.0, 2.0, 4.0, 4.5, -4.5, -7.5, 11.0, 12.5 

(8b) 

Fl : 0.5, -1.0, -2.0, -4.0, -4.5, 4.5, 7.5, -11.0, -12.5 

we recover just the leading-order diagrams listed in (7a) for the generic Dirac mass matrix 
together with the following uniquely-ordered diagrams for the complex-symmetric Majorana 
mass matrix 



M33 : 


16 3 


-1 - 126' 


- 1 


M23 : 


16 2 


-l-45 x 


-V 


M32 : 


16 3 


-1 - 126' 


- V 


M22 : 


16 2 


-l-45 x 


-V 


M 13 : 


16i 


-45 x -1' 


- 1 


M31 : 


16 3 


-1 - 126' 


- V 


M 12 : 


16i 


-45 x -1' 


- 1 


M21 : 


16 2 


-l-45 x 


-V 


Mil : 


16i 


-45 x -1' 


- 1 



±u 3 
16 2 

45 x -1 -16 2 
T26' - 1 - 16 3 



(7b) 



126 - V - 45 x -1 - 16 2 
45 z -1 - V - 45 x -16i 



Several other higher-order diagrams are allowed by the U(l)p charges given in (8a,b) and 
appear for Dll, D22, M23 and M32 with the Higgs fields ordered as follows: 



Dll : 16i -45 x -1' - 1 - 126 - 1 - 1' - 45 x -16i 

D22 : 16 2 -45 z -10'* - 1'* - 16 2 , 16 2 -1'* - 10'* - 45 z -16 2 

, (7c) 

M23 : 16 2 -45^ -V - 1 - 45 z -V - 126 - 1 - 16 3 

M32 : 16 3 -1 - T26' - V - 45 z -1 - V - 45* x -16 2 

These corrections to M23 and M32 ensure that M R is rank 3 and nonsingular, so that the 
seesaw formula can be applied. Up to this point the contributions are all complex-symmetric. 
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Additional correction terms of higher order which need not be complex-symmetric can be 
generated for the Dirac and Majorana matrix elements, if one allows additional superheavy 
fermion pairs with new U(1)f charges. Such a subset which does not destroy the pattern 
constructed above but helps to improve the numerical results given later consists of the 
following: 

F L : 1.5, -6.0, -6.5 

(8c) 

F C L : -1.5, 6.0, 6.5 
We shall enumerate the additional diagrams contributing to Dll, D12, D13, D21, D31 and 

Mil in a more detailed paper in preparation [12]. 

We assume the superheavy fermions all get massive at the same mass scale, so each 
1, 1', 45 x or 45z vertex factor can be rescaled by the same propagator mass M. As a result 
there are 14 independent parameters which can then be taken to be 

9wv u , g w v d , g w >v' d , g 12 6W u , g^e^d, 9i2&w' 

g45 x u 45x /M, g i5z u i5z /M, gi^/M, g v u v /M (9) 

045 x , 045 z , 01, 01' 

In addition, one needs the Clebsch-Gordan coefficient appearing at each vertex which can be 
read off from Table I. The algebraic contributions to each matrix element of the four Dirac 
and one Majorana matrices will be spelled out explicitly in Ref. [12]. 

One particularly good numerical choice for the parameters at the SUSY GUT scale is 
given by 

g 10 v u = 120.3, g 10 v d = 2.46, g w v' d = 0.078 GeV 

gi2%w u = 0.314, g 126 w d = -0.037, g 12& w' = 0.8 x 10 16 GeV 

g i5x U45 x /M =0.130, g 4bz u 4bz /M =0.165, g\u x jM =0.56, gyu'jM =-0.026 

045x = 35°, 45z = 01 = 01' = -5° 

(10) 

which reduces the number of independent parameters to 12 and leads to the following mass 
matrices at the SUSY GUT scale 
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M u = 



M D = 



( -0.0010 - 0.0001? 0.0053 + 0.0034i -0.0013 \ 

0.0053 + 0.0034?' 0.314 

\ -0.0013 120.3 J 

( -0.0001 -0.0104 + 0.0004?' -0.0029- 0.0015/^ 

-0.0077 + 0.0018? -0.036 0.078 

V -0.0033 - 0.0048?' 0.078 2.460 



(11a) 



(116) 



M 



N 



( 0.0030 + 0.0003?' -0.079 - 0.051? 0.0038 \ 
0.048 + 0.031? -0.942 

0.0038 120.3 j 



M E = 



V 

/ 0.0004 

0.0060 + 0.0031?' 

V -0.0009 - 0.0037?' 
/ (-.075 + .637?) x 10 9 



(11c) 



M R = 



0.0020 - O.OOlOi -0.0023 - 0.0045i \ 

0.108 0.078 (lid) 

0.078 2.460 j 

(-.139 - .117?) x 10 11 (.106 + .019i) x 10 13 \ 



(-.139- .117?) x 10 11 (.454 + .541?) x 10 12 (-.387 - .152i) x 10 



11 



lie) 



V (.106 + .019?) x 10 13 (-.387 - .152?) x 10 14 (.247 - .OUi) x 10 16 J 
in units of GeV. The Dirac mass matrix elements appear in the form ■0jC _1 M(-?/' c )i, while 
the Majorana matrix elements refer to (i/^^C -1 M(^ C ) L with ip L and (^> c )l each a member 
of one of the three families of 16's. Identical contributions also arise from the transposed 
Dirac matrices and the right-handed Majorana matrix. As such, the true Yukawa couplings 
G Y are just half the values of the gy's appearing in (9) and (10). 

The masses at the GUT scale can then be found by calculating the eigenvalues of the 
Hermitian product MM^ in each case, while the mixing matrices Vckm and V\ ev ton can 
be calculated with the projection operator technique of Jarlskog [13]. After evolving these 
quantities to the low scale, we find in the quark sector 

m u (lGeV) = 5.0 (5.1) MeV, m d (lGeV) = 7.9 (8.9) MeV 

m c (m c ) = 1.27 (1.27) GeV, m s (lGeV) = 169 (175) MeV (12a) 

m t (m t ) = 150 (165) GeV, m b (m b ) = 4.09 (4.25) GeV 



-9- 



FERMILAB-Pub-95/107-T 



where we have indicated the preferred values in parentheses. The mixing matrix is given by 

/ 0.972 0.235 0.0037 e - il24 ° \ 

-0.235 0.971 0.041 

0.999 



V, 



CRM 



V 



0.012 
-0.003i 



-0.039 
-O.OOli 



(126) 



/ 



Note that \V ub /V cb \ = 0.090 with the CP-violating phase 5 = 124°, while m d /m u = 1.59, 
m s /m d = 21.3, cf. Ref. [14]. In the lepton sector we obtain 



= 0.10 (?) x 10" 4 eV, 
= 0.29 (0.25) x 10- 2 eV, 
= 0.11 (0.10) eV, 



and 



V, 



( 0.998 
-0.036 



lept 



V 



0.043 
-0.037i 



m e = 0.44 (0.511) MeV 
m M = 99 (105.5) MeV 
m T = 1.777 (1.777) GeV 

0.050 0.038e- il22 ° \ 
0.873 0.486 
0.873 



(13a) 



-0.485 
-0.002i 



(136) 



/ 



The heavy Majorana neutrino masses are 



Mf = 0.63 x 10 9 GeV, M 2 fi = 0.39 x 10 11 GeV, = 0.25 x 10 16 GeV (13c) 

The neutrino masses and mixings are in the correct ranges to explain the nonadiabatic solar 
neutrino depletion [3] with small mixing and the atmospheric neutrino depletion [4] with 
large mixing: 



Sm 2 12 



8.5 x 10~ 6 eV 2 , 
om 2 3 = 1.2 x 10- 2 eV 2 , 



sin 2 2# 12 = 1.00 x 10~ 2 
sin 2 2# 23 =0.72 



(14) 



For our analysis, the SUSY GUT scale at which the gauge and Yukawa couplings unify 
was chosen to be A = 1.2 x 10 16 GeV. From (4a) and (11) we find that g 10 = 0.69. It is 
interesting to note that if we equate the S'O(10)-breaking and lepton number-breaking VEV, 
w', with A, we find g 12 6' = 0.67 ~ g w . Taking into account the remark following (lie), we 
note the true Yukawa couplings are G w ~ G126' — 0.33. If we further equate gi = g w ~ gi2&, 
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and Mi = A for the L^l ^-breaking VEV, we find M — 1.5 x 10 16 GeV for the masses of the 
superheavy fermions which condense with their mirrors. These values are all very reasonable. 

The 45x and 45^ VEVs appear at nearly the same scale, 2.8 x 10 15 and 3.5 x 10 15 
respectively, if one assumes the same Yukawa coupling as above. On the other hand, if these 
VEVs appear at the unification scale A the corresponding Yukawa couplings are smaller than 
those found above. In either case, a consequence of their non-orthogonal breakings is that 
577(5) is broken down to SU(3) C x SU(2) L x U(l)y at the scale in question. No further 
breaking is required until the electroweak scale and the SUSY-breaking scale are reached. 

In summary, we have constructed an 50(10) x U(l)p model of the Yukawa interactions 
with the following features: 

(i) The global U(l)p family symmetry group singles out a rather simple set of tree di- 
agrams which determines the texture of the generic Dirac and Majorana mass matrices, 
while the 50(10) group relates corresponding matrix elements of the up, down, neutrino 
and charged lepton Dirac matrices to each other. 

(ii) The dominant second and third family Yukawa interactions are renormalizable and 
arise through couplings with Higgs in the 10, 10' and 126 representations of 5O(10). The 
remaining Yukawa interactions are of higher order and require couplings of Higgs in the 
126', 1, 1', 45 x and 45 z representations which acquire VEVs near the SUSY GUT scale. 

(iii) The Higgs which acquire high scale VEVs break the 5O(10) x U(1) F symmetry down 
to the SU(3) C x SU(2) L x U(l)y standard model symmetry. 

(iv) Although this non-minimal supersymmetric model involves several Higgs represen- 
tations, the runnings of the Yukawa couplings from the GUT scale to the low-energy SUSY- 
breaking scale are controlled mainly by the contributions from the 10, as in the minimal 
supersymmetric standard model. 

(v) In terms of 12 input parameters, 15 masses (including the heavy Majorana masses) 
and 8 mixing parameters emerge. The Yukawa couplings and the Higgs VEVs are numer- 
ically feasible and successfully correlate all the quark and lepton masses and mixings in 
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the scenario which incorporates the nonadiabatic solar neutrino and atmospheric neutrino 
depletion effects. 

We shall elaborate further on the numerical details in a paper now in preparation [12]. 
Work is also underway to construct a superpotential for the model presented here. 

The authors gratefully acknowledge the Summer Visitor Program and warm hospitality 
of the Fermilab Theoretical Physics Department where much of this research was initiated 
and carried out. We thank Joseph Lykken for his continued interest and comments during 
the course of this work. The research was supported in part by the U. S. Department of 
Energy. 
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SU(5) 


VEV Directions 


Flipped SU(5) 


Assignments 




45y 


45 z 


Assignments 


u, d 


1 


1 


1 


d, u 


u c 


1 


- 4 


2 
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Table I. Couplings of the 45 VEVs to states in the 16. 



Fig. 1. Tree- level diagrams for the (a) renormalizable and (b) leading-order nonrenormal- 
izable contributions to the generic Dirac mass matrix and for the (c) 33 element of the 
Majorana mass matrix. 
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